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Abstract 



We argue that the production of forward jets in polarised deep inelastic scat- 
tering can be a useful tool for probing the double /n^(l/x) resummation effects 
^ ■ which control the polarised deep inelastic scattering for small values of the 

^^ ' Bjorken parameter x. We solve the corresponding integral equations generating 

'!::j- ■ the double /n^(l/x) resummation and calculate the differential structure function 



describing the forward jet production in the small x regime which can possibly 



On ' be probed by the polarised HERA measurements. We show that these structure 

■^ ' functions should exhibit the characteristic increase with decreasing x/xj, where 

Q^' xj denotes the longitudinal momentum fraction of the parent proton carried by 

{2L|' ^ J6t, and we quantify this increase. 



The idea of studying the deep inelastic ep scattering with an identified forward jet 
as a probe of the low x behaviour of QCD was proposed by Mueller |I[]. Since then 
it has been successfully applied to the case of unpolarised deep inelastic scattering 
(DIS) [0, ^. The Mueller's proposal was to study the unpolarised deep inelastic events 
[x, Q"^) containing an identified forward jet {xj, k'j), where the longitudinal momentum 
fraction xj of the proton carried by a forward jet and jet transverse momentum squared 
k'j were assumed to fulfill the conditions : 

Xj > X, (1) 

k'j ~ Q'- (2) 
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Here, as usual, Q^ = —q'^, where q is the four momentum transfer between electrons 
in the deep inelastic ep scattering, and x is the Bjorken parameter, i.e. x = Q'^/{2pq) 
with p denoting the four momentum of the proton. 

The first assumption xj 3> x implies that one probes the QCD dynamics in the 
low x/xj region (the quantities which are measured depend on the ratio x/xj). The 
second assumption kj ~ Q^ guarantees suppression of the standard leading order (LO) 
DGLAP evolution [|l], H, ^ and restricted penetration of the non-perturbative region by 
the small x resummation 0]. Let us recall that the small x behaviour of (unpolarised) 
deep inelastic scattering should be controlled by the QCD pomeron which in pertur- 
bative QCD is described by the Balitskij, Fadin, Kuraev, Lipatov (BFKL) equation 
0. This equation generates the leading ln{l/x) resummation, and it corresponds to 
the sum of ladder diagrams with the (reggeised) gluons along the chain. Unlike the 
LO DGLAP equations which correspond to ladder diagrams with ordered transverse 
momenta, the transverse momenta of the gluons along the BFKL ladder are not or- 
dered. In fact, the small x behaviour generated by the BFKL equation is linked with 
the diffusion of the transverse momenta both towards the infrared and the ultraviolet 
regions. It is this diffusion which is probed in the kinematical configurations which 
have comparable scales on both sides of the ladder, i.e. fcj ~ Q^ for the forward jet 
production. 

In this paper we would like to apply the idea of forward jet measurement to the case 
of polarised deep inelastic ep scattering in the region of small values of x which can be 
probed at the possible polarised HERA measurements 0. The pomeron contribution 
decouples from the polarised deep inelastic scattering, and its small x limit is sensitive 
to the novel effects of the double ln^{l/x) resummation 0, §. 

The relevant dynamical quantity in the case of jet production in polarised deep 
inelastic scattering is the differential spin structure function ^Jq^^qp- It is linked to 
the corresponding differential cross-section in the standard way: 

^^ -^-^ «^9' ,(2-,). (3) 



dxdQ'^dxjdk'j Q^ dxjdk'j 

where ag denotes the electromagnetic coupling constant, and y describes the energy 
fraction of incoming electron carried by the interacting virtual photon. In equation 
(0) we have omitted terms proportional to 7^ = AM'^x'^/Q'^, where M denotes nucleon 
mass, which are negligible at small x. The cross section in formula (^ corresponds, 
as usual, to the difference between the cross-sections for antiparallel and parallel spin 
orientations 0. Similarly as for the unpolarised case, restrictions (P, (0) applied to 
the polarised DIS forward jet events allow one to neglect the effects of DGLAP evolu- 
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Figure 1: An example of the ladder diagram contributing to the differential spin struc- 
ture function. 
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tion and to concentrate on the low x/xj behaviour of xj ^^ y^i . Low x/xj behaviour 
of differential spin structure function is dominated by double logarithmic ln^{l/x) 
(or rather ln^{xj/x)) contributions i.e. by those terms of the perturbative expansion 
which correspond to the powers of ln^{xj/x) at each order of the expansion, similarly 
as for the DIS spin dependent structure function gi [0, |]. In the following we will 
apply an approach to the double ln^{l/x) resummation based on unintegrated parton 
distributions I 



P, 12, 13 



The full contribution to the double ln'^{l/x) resummation comes from the ladder 
diagrams with quark and gluon exchanges along the ladder (see e. g. Fig. 1) and the 



non- ladder bremsstrahlung diagrams |T^. The latter ones are obtained from the ladder 
diagrams by adding to them soft bremsstrahlung gluons or soft quarks [0, H, [l^, and 



they generate the infrared corrections to the ladder contribution. However, it was 
shown in ref. |13[ that the effects of non-ladder contributions considered for DIS spin 
dependent structure function gi are visible only for x < 10~^. In what follows we 
assume x/xj > 10~^ and limit ourselves to the contributions from ladder diagrams 
which dominate in this region. 

The relevant region of phase space generating the double /n^(l/x) resummation 
from ladder diagrams corresponds to ordered k^^/xn, where k"^ and Xn denote respec- 
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Figure 2: Diagrammatic representation of formula (^ for the differential structure 
function a^j^^r-fp- describing the (forward) jet production in polarised deep inelastic 
scattering. 



tively the transverse momenta squared and longitudinal momentum fractions of the 
proton carried by partons exchanged along the ladder [^. In consequence, the ef- 
fects of the double /n^(l/x) resummation should be still visible in the kinematical 
configuration having comparable scales on both sides of the ladder (i.e. fcj ~ Q^ for 
jet production) in contrast to the LO DGLAP evolution alone which corresponds to 
ordered transverse momenta. 

The formula for the differential structure function can be written in the following 
form (see Fig. 2) : 



xj 



dxjdk'] 



X 



a 



ik]) Y: Ap,{xj, fc})AP™(0)<l>w(— , kl Q'), 

xj 



(4) 
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where Api{xj, /cj) are the (integrated) spin dependent parton distributions in the pro- 
ton. The quantities APjm(O) = ^Pim{z = 0), where APim{z) (for quarks APim{z) = 
SimAPgg) denote the LO splitting functions describing evolution of spin dependent par- 
ton densities. The indices i and m numerate quarks, antiquarks and gluons and the 
index I quarks and antiquarks respectively. Finally, ag = as/^ir. Equation (^) was 
derived, assuming strong ordering of transverse momenta {kf << k'j r^ k"^) and of the 
longitudinal momentum fractions {xm << a^j ~ xj) at the jet production vertex (see 



Fig. 2). 

The functions $m«(^5 ^j? Q^) are related to the unintegrated quark and antiquark dis- 
tributions fmi{^i k'j, kV) in the parton m, where k'j denotes the transverse momentum 
squared of the quark (antiquark): 

1 2 ^dk) r^jdx' ^x ,, ,,^ 



where 



$w(-,fcj,g^) = -efy -^j^ —u,i-,kikj), (5) 



,2 I'^J 



W' = g^(^^-lj. (7) 

We show later that the integration limits in eq. (^ are further restricted by the 
ordering of k'^/xn- In particular this ordering will provide the infrared cut-off for kj. 

It is convenient to introduce the non-singlet and singlet combination of parton 
distributions 



Nf 

mh 



^^' ~lr{W) ^'^' 



Nf 

Jm ~ 2^ Jmh WJ 

1=1 

where (e^) = -^ Hi=i ^h ^^^ ^f denotes the number of active flavours {Nj = 3). 

Restricting ourselves to ladder diagrams in the double /n^(l/^) approximation, we get 
the following equations for the functions f^^i^, kj, kj) and f^X^, k'j, kj-): 

a.(M^)AF,,(o) f ^ fj'^ ^CH^\ ki e)Q{e - ejO. m 

/m(?; kji kf) = /„j (^, kj, kf) + 



/m(?; kj, kf) — frn (C; kj, kj:) + 

(12) 



pi J ft fk^/z 71.2 

a.(M^) j yI ^ (AP,,(0)/^(e', kl e) + AP,MfL{(,', kj, e)) Q{e - k]a 



where z = ^/C,'- The inhomogeneous terms are given by the following formulae: 

/r°(e, A:}, fcj) = 5(fcJ-A:})<5(e-l)(^-l)(<5„,,„+<5^,-J, (13) 

f^'i^,klkj) = 6{kj-ej)6{^-l){6^,^ + 6^,J, (14) 

fl\^,k],kj) = 5{k]-k])5{i-l)5^,. (15) 

We have solved equations (p^ , (|ll|) and ([l^) for two choices of the scale /i^ which 
allow analytical solution of these equations: ^^ = (/cj + Q'^)/2 and /i^ = k'j/C,- In 
the first case the coupling as(/i^) does not change along the ladder, and it is a fixed 
parameter. In what follows we shall call this choice "the fixed coupling case". Let us 
remind that this choice of the scale follows the convention usually adopted in the case 
of forward jet production in the unpolarised deep inelastic scattering 0. In the second 
case /x^ = kj/C, the coupling changes along the ladder, and so we shall call this choice 
the "running coupling case". One expects that the scale //^ = kj/^ is presumably too 
large, and that this choice will lead to an underestimate of the effect. The more natural 
choice /x^ = kj does not, however, allow analytical solution. One can expect that the 
solution of equations (0) - ([I^ ) with /i^ = kj should lie between the two solutions 
corresponding to /i^ = {k'j + Q^)/2 and /i^ = kj/^. 



For the case of fixed a^ the solutions of equations ([T0|), ([TT|) and (p!2| ) for the functions 
/m^(^, fcj, kj) and /^(^, Pj, k)) take the following form : 

tm (s) kj^ kjr) = f^ (4, kj, kj) 
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+ ^('5n^,™ + 5„^,J (^ - Ij Io[2^asXyit + y)Mk'' - k'JO, (16) 



tm\iik,j,k^) — Jm{^jkj,kj) 



+ -^T=^^'^M /0[2V'«.A+ yit + y)]eikj - k]0 



asX'o 



+ -^—^A,,io)io[2^aAoy{t + y)Mkj-kJO, (17) 



where Io{z) denotes the modified Bessel function. 
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Figure 3: The differential spin structure function xjq^^^ for tlie fixed cxs case plotted 
as the function of the longitudinal momentum fraction xj carried by a jet. We show 
our predictions for three values of the transverse momentum squared fcj of the jet i.e. 
for 5 GeV^, 10 GeV^ and 20 GeV^. The calculations were performed for Q"^ = lOGeV'^ 
and X = 10~^. 



A([, Xq are the eigenvalues of matrix AP(0) 



AP,,(0) AP,,(0) 
AP,,(0) AP,,(0) 



Ao" 



AP,,(0) + AP,,(0) ± J(AP,,(0) - AP,,(0))2 + 4AP,,(0)AP,,(0) 



(21) 



and coefficients c^, c„ are defined as : 



^mg y^rnqm ' ^' 



A+-AP,,(0) 



~Omg + \Omqm + ^mqm) 



™''"^ APg3(0) ' 

Ao--AP,,(0) 



AP,,(0) 



This form is typical for the double logarithmic asymptotics (cf. 
therein). It has also been obtained for fragmentation functions [J 
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The solution for the case of running a^ = as{k'^f/C) fo^^ the non-singlet and the 
singlet components read : 

a,(fc^/OA 
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^±^]e(fc2_fc2^), (22) 
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Figure 4: The differential spin structure function Xj q^ ^^^-2 for the running as case 
plotted as the function of the longitudinal momentum fraction xj carried by a jet. 
We show our predictions for three values of the transverse momentum squared k'j of 
the jet i.e. for 5 GeV^, 10 GeV^ and 20 GeV^. The calculations were performed for 
Q2 = lOGeV^ and x = 10"^ 
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(23) 
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Having calculated the unintegrated parton distributions for both the fixed and the 
running ots case, one may express the double logarithmic contribution X)/ ^m i in (^ 
as : 
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Figure 5: Born approximation for the differential spin structure function Xj ^^ ^i^-^ 
plotted as the function of the longitudinal momentum fraction xj carried by a jet. 
We show our predictions for three values of the transverse momentum squared fcj of 
the jet i.e. for 5 GeV^, 10 GeV^ and 20 GeV^. The calculations were performed for 
Q2 = lOGeV^ and x = 10"^ 



where the integration limits x^p, Xdown, ^^p and A;^o^„ which read : 



X. 
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up 



X 



Xj, 



up 



Q\- 
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Xdo 



down 



■r' h^ 



take into account restrictions implied by the integration limits in 
functions in equations (|TB|,|T^ etc. i. e. : 



(28) 
and by the B 



< W 



X < X < Xj, 
/ 2 7 2 



(29) 
(30) 
(31) 



The latter restriction is the most important one: it follows from the ordering of the ratio 
k'^/xn along the ladder which, in turn, determines the infrared limit in the integration 
over dkj. 

We have performed numerical analysis of the low x/xj behaviour of differential 

. It was calculated by convoluting the double logarithmic 



structure function 



dx jdk' 



contribution X)/ ^mi for the fixed and running a^, respectively, with the spin dependent 
parton distributions Api{xj,kj). The parton distributions Api{xj,k'j) were obtained 
from the input distributions at the initial scale equal 1 GeV^ parametrised as in ref. 
jl^ , [13| which were evolved using the leading order Altarelli-Parisi evolution for the 
spin dependent parton densities |T^. Our results are presented in Figs. 3 and 4 for the 



fixed and running a^ case respectively. In both cases the solutions depend strongly on 
the value of the jet transverse momentum kj (cf. equations (p!6D ,([T7| ) ,(|2^ ) ,(^)). One 
may notice that the absolute value of the differential structure function exhibits very 
strong increase with the increasing value of xj (i.e. decreasing x/xj). This increase 
is a direct consequence of the double /n^(l/x) resummation summarised in equations 
(P!^),(P^),(P2D,(P^. The absolute magnitude of the differential structure function which 
corresponds to the running coupling case (i.e. to the scale //^ = kj/^ in equations 
(IIO)'(0) ^^^ dUD) is significantly smaller than that one which corresponds to the 
fixed coupling , i.e. //^ = {k'j + Q'^)/2. The former scale is presumably too large, and 
the latter case may be regarded as the more realistic one. 

We have also calculated the differential structure function in Born approximation 
which corresponds to the replacement of functions f^^i^, k'j, k"^) and /^('C) ^j? ^^) by 
inhomogeneous terms /^'^°(^, fcj, fc^) and f^{^,k'j,k'^) defined by equations (|13[) and 



(|l^). Results of this calculations are summarised in Fig. 5. We may see that the 
absolute magnitude of the differential structure function with the effects of the double 
/n^(l/x) resummation (see Figs. 3 and 4) taken into account is significantly larger 
than the absolute magnitude of the differential structure function calculated in Born 
approximation. Let us also note that the xj dependence is completely different in those 
two cases i.e. the effects of the double /n^(l/x) resummation generate increase of the 
absolute magnitude of the structure functions with increasing xj (for fixed x), while 
the absolute magnitude of the structure function calculated in Born approximation 
decreases with increasing xj, just following the xj dependence of the spin dependent 
parton distributions Api{xj,k'j). Comparison with Born approximation leads to the 
conclusion that the double logarithmic resummation effects should be clearly visible 



-2 /^ _ in-4N 



for xj > 10"^ (x = 10 

To sum up, we have estimated possible effects of the double lin?{l/x) resummation 
in the forward jet production in the polarised deep inelastic scattering at the small 
X regime which will possibly be probed at the polarised HERA. We have shown that 
the (absolute) magnitude of the differential spin dependent structure function strongly 
increases with decreasing x/xj, and that it is significantly larger than the "background" 
which corresponds to Born term with the double ln'^{l/x) resummation neglected. 
Estimate of the expected rate of the forward jet production in polarised deep inelastic 
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ep scattering taking into account the HERA acceptance is in progress. 
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